Measurements of 
Inclusive b-quark Production at 7 TeV 
with the CMS Experiment 



Abstract 



P. Bellan 

,_l ON BEHALF OF 

the CMS Collaboration 

O 

<N 

^ Padova University and INFN - 35131 Padova, Italy 

<D 

On 
D 

Oh 
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inclusive b-quark production in proton-proton collisions at y/s = 7 TeV 
are presented. The measurements are based on different methods, such 
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1 Introduction 



The study of heavy-quark production in high-energy hadronic interactions plays a 
key role in testing next-to-leading order (NLO) Quantum Chromodynamics (QCD) 
calculations. In the past, discrepancies were observed between experimental data and 
theoretical predictions, e.g. at Tevatron [Tj |2j [3J, H] and HERA [5j El [7J [8] . Substan- 
tial progress has been achieved in the understanding of heavy-quark production at 
Tevatron energies [9J, but large theoretical uncertainties still remain, mainly due to 
the dependence of the calculations on the renormalization and factorization scales. 
The observed large scale dependence of the NLO calculations is considered to be 
a symptom of large contributions from higher orders: small-x effects plj ITT] , where 
x ~ mb/v^, are possibly relevant in the low transverse momentum (pp) domain, while 
multiple-gluon radiation leads to large logarithms of / and may be important at 
high pt P2] • The resummed logarithms of £>t / ™>b at next-to-leading-logarithmic accu- 
racy have been matched to the fixed-order NLO calculation for massive quarks [TBI . 
At the non-perturbative level, the b-hadron spectrum depends strongly on the 
parametrization of the fragmentation function [13] . The b-quark production cross 
section has also been studied in the general-mass variable-flavor-number scheme [To] 
and the factorization QCD approach [TBI ITT] . 

Measurements of b-hadron production at higher energies than before, provided by 
the Large Hadron Collider (LHC), represent an important test of the new theoretical 
calculations just mentioned. Measurements of inclusive b-quark production cross 
section require identification of inclusive events in which a b-quark has been produced 
in pp collisions. In results reported here, the discrimination of the heavy quark events 
has been achieved either exclusively, reconstructing the whole decay channel of a B 
meson, or inclusively, considering hard jets. In this context, two different tagging 
techniques have been applied: reconstruction of a displaced secondary vertex, and 
analysis of the transverse momentum spectrum of an energetic muon with respect to 
the closest jet. 

2 The CMS detector 

A detailed description of the CMS detector can be found elsewhere [JS]. Some of the 
most relevant features for heavy flavor physics are summarized here. 

The core CMS apparatus is a superconducting solenoid, of 6 m internal diameter, 
providing a magnetic field of 3.8 T. Within the field volume there are the silicon pixel 
and strip tracker, the crystal electromagnetic calorimeter and the brass/scintillator 
hadron calorimeter. Muons are detected by three types of gas-ionization detectors em- 
bedded in the steel return yoke: Drift Tubes (DT), Cathode Strip Chambers (CSC), 
and Resistive Plate Chambers (RPC). The muon detectors cover a pseudorapidity 
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window < 2.4, where rj = — ln[tan(6 l /2)], where the polar angle 9 is measured 
from the z-axis, which points along the counterclockwise beam direction. The sili- 
con tracker is composed of pixel detectors (three barrel layers and two forward disks 
on each side of the detector, made of 66 million 100 x 150 /im 2 pixels) followed by 
microstrip detectors (ten barrel layers plus three inner disks and nine forward disks 
on each side of the detector, with 10 million strips of pitch between 80 and 184 /mi). 
Thanks to the strong magnetic field and the high granularity of the silicon tracker, 
the transverse momentum, p^, of the muons matched to reconstructed tracks is mea- 
sured with a resolution of about 1 % for the typical muons used in this analysis. The 
silicon tracker also provides the primary vertex position, with ~ 20 /im accuracy. 

The first level (LI) of the CMS trigger system, composed of custom hardware 
processors, uses information from the calorimeters and muon detectors to select the 
most interesting events. The High Level Trigger (HLT) further decreases the rate 
before data storage. 

3 Open beauty 

b-jets play a key role in searches for new physics beyond the Standard Model. It 
took a while to fully establish a consistency between the Tevatron data and pQCD 
predictions for b-jet production cross section. Generally speaking, b-jets cross section 
measurements are highly non-trivial, and sizable uncertainties affect both theory and 
experiment: on one side, one has to deal with a typical multi-scale problem, in which 
the center-of-mass collision energy, mass of the b-quark and the factorization and 
re-normalization scales are entangled in a subtle way; on the other hand, excellent 
performance of the tracking is required, challenging the detector full potential. 

The CMS Collaboration has published results, reviewed in the next sections, con- 
cerning two complementary measurements of b-jet production cross section. These 
results have made use of two different b-tagging techniques, were performed over 85 
nb _1 and 60 nb _1 of 2010 data, respectively [HIED], and were obtained using track or 
particle-flow jets. Typical values achieved by the CMS detector of the jet resolution 
and the energy scale uncertainty, at the time of these results were reported, were 
about 10 — 15%, and below 3%, respectively. 

3.1 b— jets with muons 

The measurement of the integrated and differential cross section of the reaction pp — > 
b + X — > fi + X has been performed on jets coming from b-quarks. The production 
of a b-quark decaying semi-leptonically is deduced by the identification of a rather 
energetic muon inside a jet, where the transverse momentum relative to the jet axis is 
quite sizable. For a muon from a b-decay, the transverse momentum relative to the jet 
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axis is on average larger than when the muon comes from light quarks; through this 
property it is hence possible to discriminate events in which b-quarks were produced. 
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Figure 1: Differential cross section (left) ~^{pp ~^ b + X — > ji + X', 1 77^ | < 2.1), and 
(right) j^ipp — » +X — > n + X'jPy > 6 GeV). The two possible muon charges are not 
distinguished and the process pp — > b + X — > jj, + X' is included. The black points are 
the CMS measurements. Vertical error bars showing the statistical error are smaller 
than the point size in most bins, the horizontal bars indicate the bin width. The 
yellow band shows the quadratic sum of statistical and systematic uncertainties. The 
systematic uncertainty (11%) of the luminosity measurement is not included. The 
solid blue line shows the MC@NLO result and the dashed blue lines illustrate the 
theoretical uncertainty as described in the text. The solid red line with dots shows 
the PYTHIA result. 



A binned log-likelihood fit is performed on the spectrum of such a quantity, called 
"pJp ei ", using template distributions provided by the simulation for b and c quarks, 
and derived from the data for gluons and lighter quarks. This latter distribution is 
dominated by hadrons misidentified as muons (mainly decay-in-flight), so they are 
reweighted by the misidentification rate measured in the data. Considering that the 
fit is not able to distinguish between light quark, gluon and charm components, these 
are merged together. The b-jet tag efficiency achieved with this technique is about 
74% at — 6 GeV, and close to 100% above 20 GeV, whereas the contamination is 
~ 7% in lowest pt bin, asymptotically decreasing towards 2% at high p^. 

The result of the inclusive production cross section for b quarks decaying into 
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muons within the kinematic range p T > 6 GeV and 1 77^ | < 2.1 is: 

a = 1.32 ± O.Ol(stat) ± 0.30(syst) ± 0.15(lumi) /ib 

where the first uncertainty is statistical, the second is systematic, and the third is 
associated with the estimation of the integrated luminosity. Fig.[T]shows the measured 
value of the cross section for the reaction pp — > b + X — > fi + X, as a function of 
the muon's pt and y. The result includes efficiency of trigger (88 ± 5)%, of muon 
reconstruction (94 ± 3)% and muon-jet association (77 ± 8)%. 

The fit stability was tested against variation of the binning, repeating the fits 
on many simulated pseudo-experiments, cross-checking the results using jets from 
particle flow and performing the fits on the impact parameter distribution. 

Table 1: Summary of systematic cross section uncertainties. The systematic uncer- 
tainty can vary depending on the muon transverse momentum and pseudorapidity as 
indicated by the range. 



source 


cross section uncertainty (%) 


Trigger efficiency 


5 


Muon reconstruction efficiency 


3 


Hadron tracking efficiency 


2 


b p™ 1 shape uncertainty 


< 21 


Background p™ 1 shape uncertainty 


2-14 


Background composition 


3-6 


Production mechanism 


2-5 


Fragmentation 


1-4 


Decay 


3 


Underlying event 


10 


Luminosity 


11 



Table [T] summarizes the systematic uncertainties entailed in this measurement. 
The uncertainties dominating the measurement are those coming from the approxi- 
mate knowledge of the signal and the background p^ 1 shape. 

Fig. [2] shows the ratio of the measured differential production cross sections in 
Pt and 77 over the predictions of the MC@NLO computation, as well as those provided 
by the Pythia [21], CASCADE MC [22] and the FONLL [23] calculations. 

3.2 b— tagging with secondary vertices 

The identification of jets coming from the hadronization of b-quarks is possible also 
through the reconstruction of secondary vertices (SV). Once a displaced SV is re- 
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Figure 2: Differential cross section (left) ^r(pp ->• b + X — > n + X' ,\r}V\ < 2.1), and 

(right) j^{pp -> b + X -> + X',^ > 6 GeV), divided by their MC@NLO predictions. 
The two possible muon charges are not distinguished; the cross section includes the process 
pp — > b + X — > fi + X. The black points are the CMS measurements. Vertical error bars 
showing the statistical error are smaller than the point size in most bins, the horizontal 
bars indicate the bin width. The yellow band shows the quadratic sum of statistical and 
systematic uncertainties. The systematic uncertainty (11%) of the luminosity measurement 
is not included. Superimposed are the FONLL result (black dotted line) with uncertainties 
(hatched band), the CASCADE result (green, dashed line) and the pythia result (red line 
with markers), divided by the MC@NLO cross section. 

constructed, different discriminators can be used to tag the jet as originating from a 
b-quark. In the analysis presented here, the discriminator adopted is a monotonic 
function of the 3D decay length. The decay length significance cut is chosen so that 
the corresponding tagging efficiency is about 60% at pp* = 100 GeV, with a con- 
tamination of ~ 0.1%. The b-tagging efficiency and the mistag rates from c or light 
jets are evaluated from simulated events and constrained by a data/MC scale factor. 
The b-tagging efficiency with the selections used in this analysis is between 6% and 
60% at pt> 18 GeV and \y\ < 2.0. The efficiency rises at higher as the b-hadron 
proper-time increases. In order to evaluate the purity of the selected sample, a fit to 
the SV mass distribution is performed, taking the shapes from simulated events, and 
letting free the relative normalizations for c and b jets, with the (small) contribution 
from light quarks fixed to the Monte Carlo expectations ("template fit"). The effi- 
ciencies estimated from MC, and the estimates of b-tagged sample purity resulting 
from fits to secondary-vertex mass from data are shown in Fig. [3] 
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Figure 3: b-tagging efficiency in different rapidity bins, as estimated on simuiated events 
(left). The b-tagged sample purity obtained using fits to secondary vertex mass (right) 

The data sample was collected using hadronic triggers with different thresholds 
of the jet pt. To merge them, the individual p^ spectra of jets have been normalized 
with the luminosity of their data taking periods, and then combined into a single 
spectrum with the jet p? bins corresponding to intervals where the triggers were fully 
efficient. The overall transverse jet energy range goes from 18 to 300 GeV, and the 
measurements have been performed in four r\ intervals. The jet energy corrections 
applied for rapidity dependence, and those for absolute scale and px dependence, 
come from real data and simulated events, respectively. 

The leading systematic uncertainties affecting the measurements are: 

• the jet energy scale of b-jets relatively to the inclusive ones (45%); 

• data-driven constraints on b-tagging efficiency (20%); 

• mistag rate for charm (34 %) and for light jets (1 - 10%). 

Fig. [4] shows the results for the production cross section measurements for b- 
jets as a function of the b-jet transverse momentum, compared with the MC@NLO 
predictions, and the ratio with the inclusive jets cross-section. While the agreement 
with Pythia and MC@NLO is reasonable, significant differences in shape are evident, 
the simulations predicting more b-jets at high p^ than what is observed. 
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Figure 4: Production cross section measurement for b-jet as a function of the b-jet 
transverse momentum compared with the MC@NLO predictions (left) and the ratio of the 
b-jet cross section over the inclusive jets (right). 

4 Conclusions 

Recent results published by the CMS Collaboration on open beauty production have 
been summarized, all of them obtained analyzing the pp collision data at 7 TeV col- 
lected in the year 2010. The studies made with two different techniques for b-tagging 
have been shown, together with the comparisons of the production cross section mea- 
sured in the data with the predictions from the available theoretical models. 
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